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Abstract 

, A summary is given of the current status and plans for gravitational-wave 

Q> ■ searches at all plausible wavelengths, from the size of the observable universe 

0^ , to a few kilometers. The anticipated scientific payoff from these searches is de- 

scribed, including expectations for detailed studies of black holes and neutron 
^ i stars, high-accuracy tests of general relativity, and hopes for the discovery of 

exotic new kinds of objects. 

in 

1 Introduction 

. There is an enormous difference between gravitational waves, and the electromag- 
■t::J- ' netic waves on which our present knowledge of the Universe is based: 

o 

' • Electromagnetic waves are oscillations of the electromagnetic field that prop- 

agate through spacetime; gravitational waves are oscillations of the "fabric" 



o 

!>■ , 

Q\ , of spacetime itself. 



Astronomical electromagnetic waves are almost always incoherent superpo- 
sitions of emission from individual electrons, atoms, or molecules. Cosmic 
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bX). gravitational waves are produced by coherent, bulk motions of huge amounts 

of mass-energy — either material mass, or the energy of vibrating, nonlinear 
^ ■ spacetime curvature. 

■ • Since the wavelengths of electromagnetic waves are small compared to their 

sources (gas clouds, stellar atmospheres, accretion disks, ...), from the waves 
we can make pictures of the sources. The wavelengths of cosmic gravitational 
waves are comparable to or larger than their coherent, bulk-moving sources, 
so we cannot make pictures from them. Instead, the gravitational waves 
are like sound; they carry, in two independent waveforms, a stereophonic, 
symphony-like description of their sources. 

• Electromagnetic waves are easily absorbed, scattered, and dispersed by mat- 
ter. Gravitational waves travel nearly unscathed through all forms and 
amounts of intervening matter |||, Q . 

• Astronomical electromagnetic waves have frequencies that begin at / ~ 10^ 
Hz and extend on upward by roughly 20 orders of magnitude. Astronomi- 
cal gravitational waves should begin at ~ 10^ Hz (1000-fold lower than the 
lowest-frequency astronomical electromagnetic waves), and should extend on 
downward from there by roughly 20 orders of magnitude. 



These enormous differences make it likely that: 

• The information brought to us by gravitational waves will be very different 
from (almost "orthogonal to") that carried by electromagnetic waves; gravi- 
tational waves will show us details of the bulk motion of dense concentrations 
of energy, whereas electromagnetic waves show us the thermodynamic state 
of optically thin concentrations of matter. 

• Most (but not all) gravitational-wave sources that our instruments detect 
will not be seen electromagnetically, and conversely, most objects observed 
electromagnetically will never be seen gravitationally. Typical electromag- 
netic sources are stellar atmospheres, accretion disks, and clouds of inter- 
stellar gas — none of which emit significant gravitational waves, while typical 
gravitational- wave sources are black holes and other exotic objects which emit 
no electromagnetic wave at all, and the cores of supernovae which are hidden 
from electromagnetic view by dense layers of surrounding stellar gas. 

• Gravitational waves may bring us great surprises. In the past, when a radi- 
cally new window has been opened onto the Universe, the resulting surprises 
have had a profound, indeed revolutionary, impact. For example, the radio 
universe, as discovered in the 1940s, 50s and 60s, turned out to be far more vi- 
olent than the optical universe; radio waves brought us quasars, pulsars, and 
the cosmic microwave radiation, and with them our first direct observational 
evidence for black holes, neutron stars, and the heat of the big bang [^. It is 
reasonable to hope that gravitational waves will bring a similar "revolution" . 

Gravitational-wave detectors and detection techniques have been under vig- 
orous development since about 1960, when Joseph Weber Q pioneered the field. 
These efforts have led to promising sensitivities in four frequency bands: 

• The Extremely Low Frequency Band (ELF), 10~^^ to 10~^^ Hz, in which the 
measured anisotropy of the cosmic microwave background radiation places 
strong limits on gravitational wave strengths — and may, in fact, have de- 
tected waves 1^ ^. The only waves expected in this band are relics of the 
big bang. 

• The Very Low Frequency Band (VLF), 10~^ to 10~^ Hz, in which Joseph 
Taylor and others have achieved remarkable gravity-wave sensitivities by the 
timing of millisecond pulsars 0|. The only expected strong sources in this 
band are processes in the very early universe — the big bang, phase tran- 
sitions of the vacuum states of quantum fields, and vibrating or colliding 
defects in the structure of spacetime, such as monopoles, cosmic strings, do- 
main walls, textures, and combinations thereof |^, ^ p!o[| . 

• The Low- Frequency Band (LF), 10~^ to 1 Hz, in which will operate the Laser 
Interferometer Space Antenna, LISA; see Sec. ^ below. This is the band of 
massive black holes (M ~ 1000 to 10^ M©) in the distant universe, and of 
other hypothetical massive exotic objects (naked singularities, soliton stars), 
as well of as binary stars (ordinary, white dwarf, neutron star, and black hole) 




Figure 1: The lines of force associated with the two polarizations of a gravitational 
wave. (From Ref. 0.) 

in our galaxy. Early universe processes should also have produced waves at 
these frequencies, as in the ELF and VLF bands. 

• The High- Frequency Band (HE), 1 to lO^Hz, in which operate earth-based 
gravitational-wave detectors such as LIGO; see Sec. ^ below. This is the 
band of stellar-mass black holes (M ~ 1 to WOOMq and of other conceivable 
stellar-mass exotic objects (naked singularities and boson stars) in the dis- 
tant universe, as well as of supernovae, pulsars, and coalescing and colliding 
neutron stars. Early universe processes should also have produced waves at 
these frequencies, as in the ELF, VLF, and LF bands. 

In this lecture I shall focus primarily on the HE and LF bands, because these 
are the ones in which observations are likely to bring us the most new information. 

2 Ground-Based Laser Interferometers: LIGO and VIRGO 

2.1 Wave Polarizations, Waveforms, and How an Interferometer 
Works 

According to general relativity theory (which I shall assume to be correct), a grav- 
itational wave has two linear polarizations, conventionally called + (plus) and x 
(cross). Associated with each polarization there is a gravitational- wave field, /i+ or 
hx , which oscillates in time and propagates with the speed of light. Each wave field 
produces tidal forces (stretching and squeezing forces) on any object or detector 
through which it passes. If the object is small compared to the waves' wavelength 
(as is the case for ground-based interferometers and resonant mass antennas), then 
relative to the object's center, the forces have the quadrupolar patterns shown in 
Fig. m. The names "plus" and "cross" are derived from the orientations of the axes 
that characterize the force patterns Q. 

A laser interferometer gravitational wave detector ("interferometer" for short) 
consists of four masses that hang from vibration-isolated supports as shown in Fig. 
^, and the indicated optical system for monitoring the separations between the 
masses ||^, Two masses are near each other, at the corner of an "L", and one 
mass is at the end of each of the L's long arms. The arm lengths are nearly equal, 
Li ~ L2 = L. When a gravitational wave, with frequencies high compared to the 
masses' ~ 1 Hz pendulum frequency, passes through the detector, it pushes the 




masses back and forth relative to each other as though they were free from their 
suspension wires, thereby changing the arm-length difference, AL = Li — L2. That 
change is monitored by laser interferometry in such a way that the variations in the 
output of the photodiode (the interferometer's output) are directly proportional to 
AL{t). 

If the waves are coming from overhead or underfoot and the axes of the + 
polarization coincide with the arms' directions, then it is the waves' + polarization 
that drives the masses, and AL{t)/L = hj^{t). More generally, the interferometer's 
output is a linear combination of the two wave fields: 

^^ = F+h+{t)+F^hAt)^h{t) . (1) 

The coefficients -F+ and F-^ are of order unity and depend in a quadrupolar manner 
on the direction to the source and the orientation of the detector The combi- 
nation h{t) of the two /I's is called the gravitational-wave strain that acts on the 
detector; and the time evolutions of h^{t), and hx{t) are sometimes called 
waveforms. 



2.2 Wave Strengths and Interferometer Arm Lengths 

The strengths of the waves from a gravitational- wave source can be estimated using 
the "Newtonian/quadrupole" approximation to the Einstein field equations. This 
approximation says that h ~ {G/c^)Q/r, where Q is the second time derivative of 
the source's quadrupole moment, r is the distance of the source from Earth (and 
G and c are Newton's gravitation constant and the speed of light). The strongest 
sources will be highly nonspherical and thus will have Q ~ ML'^, where M is their 
mass and L their size, and correspondingly will have Q ~ 2Mv'^ ~ ^-^kfn' where v 
is their internal velocity and E^f^ is the nonspherical part of their internal kinetic 
energy. This provides us with the estimate 



i.e., h is about 4 times the gravitational potential produced at Earth by the mass- 
equivalent of the source's nonspherical, internal kinetic energy — made dimension- 
less by dividing by (? . Thus, in order to radiate strongly, the source must have a 
very large, nonspherical, internal kinetic energy. 

The best known way to achieve a huge internal kinetic energy is via gravity; 
and by energy conservation (or the virial theorem), any gravitationally-induced 
kinetic energy must be of order the source's gravitational potential energy. A huge 
potential energy, in turn, requires that the source be very compact, not much larger 
than its own gravitational radius. Thus, the strongest gravity-wave sources must be 
highly compact, dynamical concentrations of large amounts of mass (e.g., colliding 
and coalescing black holes and neutron stars). 

Such sources cannot remain highly dynamical for long; their motions will be 
stopped by energy loss to gravitational waves and/or the formation of an all- 
encompassing black hole. Thus, the strongest sources should be transient. More- 
over, they should be very rare — so rare that to see a reasonable event rate will 
require reaching out through a substantial fraction of the Universe. Thus, just as 
the strongest radio waves arriving at Earth tend to be extragalactic, so also the 
strongest gravitational waves are likely to be extragalactic. 

For highly compact, dynamical objects that radiate in the high-frequency band, 
e.g. colliding and coalescing neutron stars and stellar-mass black holes, the inter- 
nal, nonspherical kinetic energy E^^/(? is of order the mass of the Sun; and, 
correspondingly, Eq. (||) gives h ~ 10"'^^ for such sources at the Hubble distance 
(3000 Mpc, i.e., 10^° light years); h ~ 10"^^ at 200 Mpc (a best-guess distance 
for several neutron-star coalescences per year; see Section |]2D, /i ~ lO-^o at the 
Virgo cluster of galaxies (15 Mpc); and h ~ 10~^^ in the outer reaches of our 
own Milky Way galaxy (20 kpc). These numbers set the scale of sensitivities that 
ground-based interferometers seek to achieve: h ~ 10^^^ to 10^^^. 

When one examines the technology of laser interferometry, one sees good pro- 
spects to achieve measurement accuracies AL ~ 10~^^ cm (1/1000 the diameter 
of the nucleus of an atom). With such an accuracy, an interferometer must have 
an arm length L = /S.L/h ~ 1 to 10 km, in order to achieve the desired wave 
sensitivities, 10~^^ to 10^^^. This sets the scale of the interferometers that are now 
under construction. 

2.3 LIGO, VIRGO, and the International Gravitational Wave Net- 
work 

Interferometers are plagued by non-Gaussian noise, e.g. due to sudden strain re- 
leases in the wires that suspend the masses. This noise prevents a single interfer- 
ometer, by itself, from detecting with confidence short-duration gravitational- wave 
bursts (though it might be possible for a single interferometer to search for the peri- 
odic waves from known pulsars). The non-Gaussian noise can be removed by cross 
correlating two, or preferably three or more, interferometers that are networked 
together at widely separated sites. 

The technology and techniques for such interferometers have been under devel- 
opment for 25 years, and plans for km-scale interferometers have been developed 
over the past 15 years. An international network consisting of three km-scale inter- 



Figure 3: Artist's conception of one of the LIGO interferometers. [Courtesy the 
LIGO Project.] 

ferometers, at three widely separated sites, is now in the early stages of construc- 
tion. It includes two sites of the American LIGO Project ("Laser Interferometer 
Gravitational Wave Observatory") and one site of the French/Italian VIRGO 
Project (named after the Virgo cluster of galaxies) p^ . 

LIGO will consist of two vacuum facilities with 4-kilometer-long arms, one in 
Hanford, Washington (in the northwestern United States; Fig. ^) and the other 
in Livingston, Louisiana (in the southeastern United States). These facilities are 
designed to house many successive generations of interferometers without the ne- 
cessity of any major facilities upgrade; and after a planned future expansion, they 
will be able to house several interferometers at once, each with a different optical 
configuration optimized for a different type of wave (e.g., broad-band burst, or 
narrow-band periodic wave, or stochastic wave). 

The LIGO facilities and their first interferometers are being constructed by a 
team of about 80 physicists and engineers at Caltech and MIT, led by Barry Barish 
(the PI), Gary Sanders (the Project Manager), Albert Lazzarini, Rai Weiss, Stan 
Whitcomb, and Robbie Vogt (who directed the project during the pre-construction 
phase). Other research groups from many different universities are contributing 
to R&D for enhancements of the first interferometers, or are computing theo- 
retical waveforms for use in data analysis, or are developing data analysis tech- 
niques. These groups are linked together by an organization called the Ligo 
Research Community. For further details, see the LIGO World Wide Web Site, 
http:/ /www. ligo. caltech.edu/ . 

The VIRGO Project is building one vacuum facility in Pisa, Italy, with 3- 
kilometer-long arms. This facility and its first interferometers are a collaboration 
of more than a hundred physicists and engineers at the INFN (Frascati, Napoli, 
Perugia, Pisa), LAL (Orsay), LAPP (Annecy), LOA (Palaiseau), IPN (Lyon), ES- 
PCI (Paris) , and the University of Illinois (Urbana) , under the leadership of Alain 
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Figure 4: LIGO's projected broad-band noise sensitivity to bursts /isB (Refs. |Tl| , 
|13[) compared with the strengths of the waves from several hypothesized sources. 
The signal to noise ratios are \/2 higher than in Ref. [^] because of a factor 2 error 
in Eq. (29) of Ref . §. 



Brillet and Adalberto Giazotto. 

Both LIGO and VIRGO are scheduled for completion in the late 1990s, and 
their first gravitational- wave searches are likely to be performed in 2001. Figure ^ 
shows the design sensitivities for LIGO's first interferometers (which are expected to 
be achieved in 2001) |jll|] and for enhanced versions of those interferometers (which 
is expected to be operating five years or so later) [^], along with a benchmark 
sensitivity goal for a more advanced interferometer that may operate in the LIGO 
vacuum system some years later |11, |l^] . 

For each interferometer, the quantity shown is the "sensitivity to bursts" that 
come from a random direction, /isb(/) [0]- This hsB is about 5 times worse than 
the rms noise level in a bandwidth A/ ~ / for waves with a random direction 
and polarization, and about 5\/5 ~ 11 worse than the the rms noise level /irms for 
optimally directed and polarized waves. (In much of the literature, the quantity 
plotted is hrms — ^sb/11-) 

This hsB is to be compared with the "characteristic amplitude" hdf) = h^/n of 
the waves from a source; here h is the waves' amplitude when they have frequency 
/, and n is the number of cycles the waves spend in a bandwidth A/ ~ / near 
frequency / [||, |ll|. Any source with he > /isB should be detectable with high 
confidence, even if it arrives only once per year. As examples. Fig. ^ shows the 
characteristic amplitudes he of several binary systems made of IAMq neutron stars 



("NS") or 10, 25, or 30 Mq black holes ("BH"), which spiral together and collide 
under the driving force of gravitational radiation reaction. As the bodies spiral 
inward, their waves sweep upward in frequency (rightward across the figure along 
the dashed lines). From the figure we see that LIGO's first interferometers should 
be able to detect waves from the inspiral of a NS/NS binary out to a distance of 
30Mpc and from the final collision and merger of a 25M0/25-/Vf0 BH/BH binary 
out to about BOOMpc. 

LIGO alone, with its two sites which have parallel arms, will be able to detect 
an incoming gravitational wave, measure one of its two waveforms, and (from the 
time delay between the two sites) locate its source to within a ~ 1° wide annulus 
on the sky. LIGO and VIRGO together, operating as a coordinated international 
network, will be able to locate the source (via time delays plus the interferometers' 
beam patterns) to within a 2-dimensional error box with size between several tens 
of arcminutes and several degrees, depending on the source direction and on the 
amount of high-frequency structure in the waveforms. They will also be able to 
monitor both waveforms /i+(t) and hx{t) (except for frequency components above 
about IkHz and below about 10 Hz, where the interferometers' noise becomes 
severe) . 

A British/German group is constructing a 600- meter interferometer called GEO 
600 near Hanover Germany |lj], and Japanese groups, a 300- meter interferometer 
called TAMA near Tokyo. GEO600 may be a significant player in the interferomet- 
ric network in its early years (by virtue of cleverness and speed of construction) , but 
because of its short arms it cannot compete in the long run. GEO600 and TAMA 
will both be important development centers and testbeds for interferometer tech- 
niques and technology, and in due course they may give rise to kilometer-scale 
interferometers like LIGO and VIRGO, which could significantly enhance the net- 
work's all-sky coverage and ability to extract information from the waves. 



2.4 Narrow-Band Detectors 

At frequencies / ^ lOOOHz, the interferometers' photon shot noise becomes a se- 
rious obstacle to wave detection. However, narrow-band detectors specially opti- 
mized for kHz frequencies show considerable promise. These include a special "dual 
recycled interferometer" [|l5|, and huge spherical or icosahedral resonant-mass de- 
tectors [16 1 that are modern variants of Joseph Weber's original "bar" detector 



3 LISA: The Laser Interferometer Space Antenna 

The Laser Interferometer Space Antenna (LISA) is the most promising detector 
for gravitational-waves in the low-frequency band, 10~^-1 Hz (10,000 times lower 
than the LIGO/VIRGO high-frequency band). 

LISA was originally conceived (under a different name) by Peter Bender of the 
University of Golorado, and is being developed by an international team led by 
Karsten Danzmann of the University of Hanover (German) and James Hough of 
Glasgow University (UK). The European Space Agency tentatively plans to fly it 
sometime in the 2014-2018 time frame as part of ESA's Horizon 2000-1- Program 




Figure 5: LISA's orbital configuration, with LISA magnified in arm length by a 
factor ~ 10 relative to the solar system. 

of large space missions. If NASA contributes significantly to LISA, then the flight 
could be sooner. 

As presently conceived, LISA will consist of six compact, drag-free spacecraft 
(i.e. spacecraft that are shielded from buffeting by solar wind and radiation pres- 
sure, and that thus move very nearly on geodesies of spacetime). All six spacecraft 
would be launched simultaneously by a single Ariane rocket. They would be placed 
into the same heliocentric orbit as the Earth occupies, but would follow 20° behind 
the Earth; cf. Fig. The spacecraft would fly in pairs, with each pair at the 
vertex of an equilateral triangle that is inclined at an angle of 60° to the Earth's 
orbital plane. The triangle's arm length would be 5 million km (10^ times larger 
than LIGO's arms!). The six spacecraft would track each other optically, using 
one- Watt YAG laser beams. Because of diffraction losses over the 5 X lO^km arm 
length, it is not feasible to reflect the beams back and forth between mirrors as is 
done with LIGO. Instead, each spacecraft will have its own laser; and the lasers 
will be phase locked to each other, thereby achieving the same kind of phase- 
coherent out-and-back light travel as LIGO achieves with mirrors. The six-laser, 
six-spacecraft configuration thereby functions as three, partially independent but 
partially redundant, gravitational-wave interferometers. 

Figure |6| depicts the expected noise and sensitivity of LISA in the same language 
as we have used for LIGO (Fig. 1): /isB = ^V^h rms is the sensitivity for high- 
confidence detection {S/N = 5) of a broad-band burst coming from a random 
direction, assuming Gaussian noise. 

At frequencies / ^ lO^^Hz, LISA's noise is due to photon counting statistics 
(shot noise). The sensitivity curve steepens at / ~ 3 x 10~^Hz because at larger / 
than that, the waves' period is shorter than the round-trip light travel time in one of 
LISA's arms. Below lO-^Hz, the noise is due to buffeting-induced random motions 
of the spacecraft that are not being properly removed by the drag-compensation 
system. Notice that, in terms of dimensionless amplitude, LISA's sensitivity is 
roughly the same as that of LIGO's first interferometers (Fig. Q), but at 100,000 
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Figure 6: LISA's projected sensitivity to bursts hsB, compared with the strengths 
of the waves from several low- frequency sources. 

times lower frequency. Since the waves' energy flux scales as f'^h?-, this corresponds 
to 10^° better energy sensitivity than LIGO. 

LISA can detect and study, simultaneously, a wide variety of different sources 
scattered over all directions on the sky. The key to distinguishing the different 
sources is the different time evolution of their waveforms. The key to determining 
each source's direction, and confirming that it is real and not just noise, is the 
manner in which its waves' amplitude and frequency are modulated by LISA's 
complicated orbital motion — a motion in which the interferometer triangle rotates 
around its center once per year, and the interferometer plane precesses around the 
normal to the Earth's orbit once per year. Most sources will be observed for a year 
or longer, thereby making full use of these modulations. 



4 Binaries Made of Neutron Stars, Black Holes, and 
Other Exotic Objects: Inspiral, Collision, and Coa- 
lescence in the HF Band 

The best understood of all gravitational-wave sources are coalescing, compact bina- 
ries composed of neutron stars (NS) and black holes (BH). These NS/NS, NS/BH, 
and BH/BH binaries may well become the "bread and butter" of the LIGO/VIRGO 
and LISA diets. 



The Hulse- Taylor 18| binary pulsar, PSR 1913+16, is an example of a 
NS/NS binary whose waves could be measured by LIGO/VIRGO, if we were to 
wait long enough. 



At present PSR 1913+16 has an orbital frequency of about 1/(8 hours) and 
emits its waves predominantly at twice this frequency, roughly 10~^ Hz, which is 
in LISA'S low- frequency band (cf. Fig. ^). As a result of their loss of orbital energy 
to gravitational waves, the PSR 1913+16 NS's are gradually spiraling inward at 
a rate that agrees with general relativity's prediction to within the measurement 
accuracy (a fraction of a percent) |l^ — a remarkable but indirect confirmation that 
gravitational waves do exist and are correctly described by general relativity. For 
this discovery and other aspects of PSR 1913+16, Princeton's Russell Hulse and 
Joseph Taylor have been awarded the Nobel Prize. 

If we wait roughly 10^ years, this inspiral will bring the waves into the LIGO/ 
VIRGO high-frequency band. As the NS's continue their inspiral, the waves will 
then sweep upward in frequency, over a time of about 15 minutes, from 10 Hz to 
~ 10^ Hz, at which point the NS's will collide and merge. It is this last 15 minutes 
of inspiral, with ~ 16, 000 cycles of waveform oscillation, and the final merger, that 
LIGO/VIRGO seeks to monitor. 



4.1 Wave Strengths Compared to LIGO Sensitivities 

Figure Q above compares the projected sensitivities of interferometers in LIGO 
[ |11| with the wave strengths from the last few minutes of inspiral of BH/BH, 
NS/BH, and NS/NS binaries at various distances from Earth. Notice that the 
signal strengths in Fig. |^ are in good accord with our rough estimates based on 
Eq. (^); at the endpoint (right end) of each inspiral, the number of cycles n spent 
near that frequency is of order unity, so the quantity plotted, he ^ h^/n, is about 
equal to h — and for a NS/NS binary at 200 Mpc is roughly 10~^^, as we estimated 



in Section 2.2 



4.2 Coalescence Rates 

Such final coalescences are few and far between in our own galaxy: about one every 
100,000 years, according to 1991 estimates by Phinney ||l^ and by Narayan, Piran, 
and Shemi [20|, based on the statistics of binary pulsar searches in our galaxy which 



found three that will coalesce in less than lO^'^ years. Extrapolating out through the 
universe, Phinney and Narayan et. al. infer that to see several NS/NS coalescences 
per year, LIGO/VIRGO will have to look out to a distance of about 200 Mpc (give 
or take a factor ~ 2); cf. Fig. ^. Since these estimates were made, the binary pulsar 
searches have been extended through a significantly larger volume of the galaxy 
than before, and no new ones with coalescence times ^ lO^'' years have been found. 
This would drive the estimated event rate downward, except that revisions of other 
aspects of the Phinney /Narayan analyses have compensated, leaving the binary- 
pulsar-search-based best estimates unchanged, at several events per year out to 200 

Mpc m]. 

A rate of one every 100,000 years in our galaxy is ~ 100 times smaller than the 
birth rate of the NS/NS binaries' progenitors: massive, compact, main-sequence 
binaries |2^. Therefore, either 99 per cent of progenitors fail to make it to 
the NS/NS state (e.g., because of binary disruption during a supernova or forming 
TZO's), or else they do make it, but they wind up as a class of NS/NS binaries 



that has not yet been discovered in any of the pulsar searches. Several experts on 



binary evolution have argued for the latter |2^, 24 1: most NS/NS binaries, they 
suggest, may form with such short orbital periods that their lifetimes to coalescence 
are significantly shorter than normal pulsar lifetimes (~ 10^ years); and with such 
short lifetimes, they have been missed in pulsar searches. By modeling the evolution 
of the galaxy's binary star population, these binary experts arrive at best estimates 
as high as 3 X 10~^ coalescences per year in our galaxy, corresponding to several per 



year out to 60 Mpc distance |22], though more recent and more conservative models 



[25 1 give results more nearly in accord with the binary pulsar searches, several per 



year out to 200Mpc. Phinney ll9[| describes other plausible populations of NS/NS 
binaries that could increase the event rate, and he argues for "ultraconservative" 
lower and upper limits of 23 Mpc and lOOOMpc for how far one must look to see 
several coalescence per year. 

By comparing these rate estimates with the signal strengths in Fig. ^, we see 
that: (i) The first interferometers in LIGO/VIRGO (ca. 2001) have a possibility 
but not high probability of seeing NS/NS coalescences, (ii) Enhanced interfer- 
ometers (mid 2000's) can be fairly confident of seeing them; the conservatively 
estimated event rate is ~ 3 x (300/200)'^ ~ 10/year. (iii) Advanced interferometers 
are almost certain of seeing them (the requirement that this be so was one factor 
that forced the LIGO/VIRGO arm lengths to be so long, several kilometers). 

We have no good observational handle on the coalescence rate of NS/BH or 
BH/BH binaries. However, theory suggests that their progenitors might not disrupt 
during the stellar collapses that produce the NS's and BH's, so their coalescence 
rate could be about the same as the birth rate for their progenitors: ~ 1/100, 000 
years in our galaxy. This suggests that within 200 Mpc distance there might be 
several NS/BH or BH/BH coalescences per year. [H, |0|, m, ||. This estimate 
should be regarded as a plausible upper limit on the event rate and lower limit on 
the distance to look |l^, pc|] . 

If this estimate is correct, then NS/BH and BH/BH binaries will be seen before 
NS/NS, and might be seen by the first LIGO/VIRGO interferometers or soon 
thereafter |2^; cf. Fig. ^. However, this estimate is far less certain than the (rather 
uncertain) NS/NS estimates! 

Once coalescence waves have been discovered, each further improvement of 
sensitivity by a factor 2 will increase the event rate by 2^ ~ 10. Assuming a rate of 
several NS/NS per year at 200 Mpc, the advanced interferometers of Fig. ^should 
see ~ 100 per year. 

4.3 Inspiral Waveforms and the Information They Can Bring 

Neutron stars and black holes have such intense self gravity that it is exceedingly 
difficult to deform them. Correspondingly, as they spiral inward in a compact 
binary, they do not gravitationally deform each other significantly until several 



orbits before their final coalescence [27,|2g]. This means that the inspiral waveforms 
are determined to high accuracy by only a few, clean parameters: the masses 
and spin angular momenta of the bodies, and the initial orbital elements (i.e. the 
elements when the waves enter the LIGO/VIRGO band). 

Though tidal deformations are negligible during inspiral, relativistic effects can 
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Figure 7: Waveforms from the inspiral of a compact binary, computed using New- 
tonian gravity for the orbital evolution and the quadrupole-moment approximation 
for the wave generation. (From Ref. [ll 



be very important. If, for the moment, we ignore the relativistic effects — i.e., if we 
approximate gravity as Newtonian and the wave generation as due to the binary's 
oscillating quadrupole moment |jl], then the shapes of the inspiral waveforms h+{t) 
and hx{t) are as shown in Fig. |7. 

The left-hand graph in Fig. ^ shows the waveform increasing in amplitude and 
sweeping upward in frequency (i.e., undergoing a "chirp") as the binary's bodies 
spiral closer and closer together. The ratio of the amplitudes of the two polariza- 
tions is determined by the inclination t of the orbit to our line of sight (lower right 
in Fig. ^). The shapes of the individual waves, i.e. the waves' harmonic content, 
are determined by the orbital eccentricity (upper right). (Binaries produced by 
normal stellar evolution should be highly circular due to past radiation reaction 
forces, but compact binaries that form by capture events, in dense star clusters that 
might reside in galactic nuclei |2£], could be quite eccentric.) If, for simplicity, the 
orbit is circular, then the rate at which the frequency sweeps or "chirps", df /dt [or 
equivalently the number of cycles spent near a given frequency, n = f'^{df/dt)~'^] 
is determined solely, in the Newtonian/quadrupole approximation, by the binary's 
so-called chirp mass, = {MiM2f/^ /{Mi + M2)^/^ (where Mi and M2 are the 
two bodies' masses). The amplitudes of the two waveforms are determined by the 
chirp mass, the distance to the source, and the orbital inclination. Thus (in the 
Newtonian/quadrupole approximation), by measuring the two amplitudes, the fre- 
quency sweep, and the harmonic content of the inspiral waves, one can determine 
as direct, resulting observables, the source's distance, chirp mass, inclination, and 
eccentricity [|^, |3^ . 

As in binary pulsar observations so also here, relativistic effects add further 
information: they influence the rate of frequency sweep and produce waveform 
modulations in ways that depend on the binary's dimensionless ratio 77 = ///M 
of reduced mass // = MiM2/(Mi + M2) to total mass M = Mi + M2 and on 
the spins of the binary's two bodies. These relativistic effects are reviewed and 



discussed at length in Refs. |32, 33|. Two deserve special mention: (i) As the waves 



emerge from the binary, some of them get backscattered one or more times off the 
binary's spacetime curvature, producing wave tails. These tails act back on the 
binary, modifying its inspiral rate in a measurable way. (ii) If the orbital plane is 
inclined to one or both of the binary's spins, then the spins drag inertial frames in 
the binary's vicinity (the "Lense-Thirring effect"), this frame dragging causes the 



orbit to precess, and the precession modulates the waveforms |35|| . 

Remarkably, the relativistic corrections to the frequency sweep — tails, spin- 
induced precession and others — will be measurable with rather high accuracy, 
even though they are typically ^ 10 per cent of the Newtonian contribution, and 
even though the typical signal to noise ratio will be only ~ 9. The reason is as 



follows H, |33, ||: 



The frequency sweep will be monitored by the method of "matched filters" ; in 
other words, the incoming, noisy signal will be cross correlated with theoretical 
templates. If the signal and the templates gradually get out of phase with each 
other by more than ~ 1/10 cycle as the waves sweep through the LIGO/ VIRGO 
band, their cross correlation will be significantly reduced. Since the total number 
of cycles spent in the LIGO/VIRGO band will be ~ 16,000 for a NS/NS binary, 
~ 3500 for NS/BH, and ~ 600 for BH/BH, this means that LIGO/VIRGO should 
be able to measure the frequency sweep to a fractional precision ^ 10^^, compared 
to which the relativistic effects are very large. (This is essentially the same method 
as Joseph Taylor and colleagues use for high-accuracy radio- wave measurements of 
relativistic effects in binary pulsars [|l8|.) 

Preliminary analyses, using the theory of optimal signal processing, predict 
the following typical accuracies for LIGO/VIRGO measurements based solely on 
the frequency sweep (i.e., ignoring modulational information) |^]: (i) The chirp 
mass Mc will typically be measured, from the Newtonian part of the frequency 
sweep, to ~ 0.04% for a NS/NS binary and 0.3% for a system containing at least 
one BH. (ii) If we are confident (e.g., on a statistical basis from measurements of 
many previous binaries) that the spins are a few percent or less of the maximum 
physically allowed, then the reduced mass fi will be measured to ~ 1% for NS/NS 
and NS/BH binaries, and ~ 3% for BH/BH binaries. (Here and below NS means 
a ~ IAMq neutron star and BH means a ~ lOM© black hole.) (iii) Because 
the frequency dependences of the (relativistic) fi effects and spin effects are not 
sufficiently different to give a clean separation between // and the spins, if we have 
no prior knowledge of the spins, then the spin//i correlation will worsen the typical 
accuracy of ^ by a large factor, to ~ 30% for NS/NS, ~ 50% for NS/BH, and a 
factor ~ 2 for BH/BH. These worsened accuracies might be improved somewhat by 
waveform modulations caused by the spin- induced precession of the orbit |34, 35 1, 
and even without modulational information, a certain combination of fi and the 
spins will be determined to a few per cent. Much additional theoretical work is 
needed to firm up the measurement accuracies. 

To take full advantage of all the information in the inspiral waveforms will 
require theoretical templates that are accurate, for given masses and spins, to 
a fraction of a cycle during the entire sweep through the LIGO/VIRGO band. 
Such templates are being computed by an international consortium of relativity 
theorists (Blanchet and Damour in France, Iyer in India, Will and Wiseman in 
the U.S.) p^l, using post-Newtonian expansions of the Einstein field equations. 



This enterprise is rather hke computing the Lamb shift to high order in powers of 
the fine structure constant, for comparison with experiment. Cutler and Flanagan 
[40 1 have estimated the order to which the computations must be carried in order 



that systematic errors in the theoretical templates will not significantly impact 
the information extracted from the LIGO/VIRGO observational data. The answer 
appears daunting: radiation-reaction effects must be computed to three full post- 
Newtonian orders [six orders in vjc =(orbital velocity)/(speed of light)] beyond 
the leading-order radiation reaction, which itself is 5 orders m. vjc beyond the 
Newtonian theory of gravity — though by judicious use of Fade approximates, these 



requirements might be relaxed [41[. 



It is only about ten years since controversies over the leading-order radiation 



reaction [42[ were resolved by a combination of theoretical techniques and binary 
pulsar observations. Nobody dreamed then that LIGO/VIRGO observations may 
require pushing post-Newtonian computations onward from 0\{vjcf \ to 0\{vjc)^^\. 
This requirement epitomizes a major change in the field of relativity research: At 
last, 80 years after Einstein formulated general relativity, experiment has become 
a major driver for theoretical analyses. 

Remarkably, the goal of 0[(z;/c)^^] is achievable. The most difficult part of 
the computation, the radiation reaction, has been evaluated to 0\(vlcf\ beyond 
Newton by the French/Indian/ American consortium [39[ and 0[(f/c)^'^] is coming 



under control. 

These high-accuracy waveforms are needed only for extracting information from 
the inspiral waves, after the waves have been discovered; they are not needed for 
the discovery itself. The discovery is best achieved using a different family of the- 
oretical waveform templates, one that covers the space of potential waveforms in a 
manner that minimizes computation time instead of a manner that ties quantita- 
tively into general relativity theory []3^ , Such templates are in the early stage 
of development. 

4.4 Testing GR, Measuring the Cosmological Universe, Mapping 
Black Holes, and Searching for Exotic Objects 

LIGO/VIRGO observations of compact binary inspiral have the potential to bring 
us far more information than just binary masses and spins: 

• They can be used for high-precision tests of general relativity. In scalar-tensor 
theories (some of which are attractive alternatives to general relativity [p4[), 
radiation reaction due to emission of scalar waves places a unique signature on 
the gravitational waves that LIGO/VIRGO would detect — a signature that 
can be searched for with high precision [|45[[ . 

• They can be used to measure the Universe's Hubble constant, deceleration 



parameter, and cosmological constant [^, |31|, |46|, |43]. The keys to such 
measurements are that: (i) Advanced interferometers in LIGO/VIRGO will 
be able to see NS/NS out to cosmological redshifts z ~ 0.3, and NS/BH 
out to z ~ 2. (ii) The direct observables that can be extracted from the 
observed waves include the source's luminosity distance tl (measured to ac- 
curacy ~ 10 per cent in a large fraction of cases), and its direction on the 



sky (to accuracy ~ 1 square degree) — accuracies good enough that only one 
or a few electromagneticahy-observed clusters of galaxies should fall within 
the 3-dimensional gravitational error boxes, thereby giving promise to joint 
gravitational/electromagnetic statistical studies, (iii) Another direct gravita- 
tional observable is (1 + z)M where z is redshift and M is any mass in the 
system (measured to the accuracies quoted above). Since the masses of NS's 
in binaries seem to cluster around 1.4Mq, measurements of (1 + z)M can 
provide a handle on the redshift, even in the absence of electromagnetic aid. 

For a NS or small BH spiraling into a massive ~ 50 to 500Mq compact central 
body, the inspiral waves will carry a "map" of the massive body's spacetime 
geometry. This map can be used to determine whether the massive body is 
a black hole (in which case the geometry must be that of Kerr, with all its 
features uniquely determined by its mass and angular momentum) or is some 
other kind of exotic compact object, e.g. a soliton star or naked singularity 



As we shall see in Section 6.2 below, this type of black-hole study 



and search for exotic objects can be carried out with much higher precision 
by LISA at large masses and low frequencies, than by LIGO/VIRGO at low 
masses and high frequencies. 

4.5 Coalescence Waveforms and their Information 

The waves from the binary's final coalescence can bring us new types of information. 

BH/BH Coalescence 

In the case of a BH/BH binary, the coalescence will excite large-amplitude, 
highly nonlinear vibrations of spacetime curvature near the coalescing black-hole 
horizons — a phenomenon of which we have very little theoretical understanding 
today. Especially fascinating will be the case of two spinning black holes whose 
spins are not aligned with each other or with the orbital angular momentum. Each 
of the three angular momentum vectors (two spins, one orbital) will drag space in 
its vicinity into a tornado-like swirling motion — the general relativistic "dragging of 
inertial frames," so the binary is rather like two tornados with orientations skewed 
to each other, embedded inside a third, larger tornado with a third orientation. 
The dynamical evolution of such a complex configuration of coalescing spacetime 
warpage (as revealed by its emitted waves) might bring us surprising new insights 
into relativistic gravity Moreover, if the sum of the BH masses is fairly large. 



~ 40 to 2OOM0 , then the waves should come off in a frequency range / ~ 40 to 200 
Hz where the LIGO /VIRGO broad-band interferometers have their best sensitivity 
and can best extract the information the waves carry. 

To get full value out of such wave observations will require ||2^ having theo- 
retical computations with which to compare them. There is no hope to perform 
such computations analytically; they can only be done as supercomputer simula- 
tions. The development of such simulations is being pursued by several research 
groups, including an eight-university American consortium of numerical relativists 
and computer scientists called the Two-Black- Hole Grand Challenge Alliance ||5^ I 
have a bet with Richard Matzner, the lead PI of this alliance, that LIGO/VIRGO 



will discover waves from such coalescences with misaligned spins before the Alliance 
is able to compute them. 



NS/NS Coalescence 

The final coalescence of NS /NS binaries should produce waves that are sensitive 
to the equation of state of nuclear matter, so such coalescences have the potential 



to teach us about the nuclear equation of state [11, 32]. In essence, LIGO/ VIRGO 
will be studying nuclear physics via the collisions of atomic nuclei that have nucleon 
numbers A ~ 10^^ — somewhat larger than physicists are normally accustomed to. 
The accelerator used to drive these nuclei up to the speed of light is the binary's 
self gravity, and the radiation by which the details of the collisions are probed is 
gravitational. 

Unfortunately, the final NS/NS coalescence will emit its gravitational waves 
in the kHz frequency band (800Hz ^ f ^ 2500IIz) where photon shot noise will 
prevent them from being studied by the standard, "workhorse," broad-band in- 
terferometers of Fig. ^. However, a specially configured ( "dual-recycled" ) interfer- 
ometer, which could have enhanced sensitivity in the kHz region at the price of 
reduced sensitivity elsewhere, may be able to measure the waves and extract their 
equation of state information, as might massive, spherical, resonant-mass detectors 
[]32| , Such measurements will be very difficult and are likely only when the 
LIGO/VIRGO network has reached a mature stage. 

A number of research groups ||5^ are engaged in numerical astrophysics simu- 
lations of NS/NS coalescence, with the goal not only to predict the emitted gravi- 
tational waveforms and their dependence on equation of state, but also (more im- 
mediately) to learn whether such coalescences might power the 7-ray bursts that 
have been a major astronomical puzzle since their discovery in the early 1970s. 

NS/NS coalescence is a popular explanation for the 7-ray bursts because (i) 
the bursts are isotropically distributed on the sky, (ii) they have a distribution of 
number versus intensity that suggests they might lie at near-cosmological distances, 
and (iii) their event rate is roughly the same as that predicted for NS/NS coales- 
cence (~ 1000 per year out to cosmological distances, if they are cosmological) . 
If LIGO/VIRGO were now in operation and observing NS/NS inspiral, it could 
report definitively whether or not the 7-bursts are produced by NS/NS binaries; 
and if the answer were yes, then the combination of 7-burst data and gravitational- 
wave data could bring valuable information that neither could bring by itself. For 
example, it would reveal when, to within a few msec, the 7-burst is emitted rel- 
ative to the moment the NS's first begin to touch; and by comparing the 7 and 
gravitational times of arrival, we might test whether gravitational waves propagate 
with the speed of light to a fractional precision of ~ O.Olsec/3 x 10^ lyr = 10~^^. 

NS/BH Coalescence 

A NS spiraling into a BH of mass M ^ lOM© should be swallowed more or less 
whole. However, if the BH is less massive than roughly IOMq, and especially if it 
is rapidly rotating, then the NS will tidally disrupt before being swallowed. Little 
is known about the disruption and accompanying waveforms. To model them with 
any reliability will likely require full numerical relativity, since the circumferences 



of the BH and NS will be comparable and their physical separation at the moment 
of disruption will be of order their separation. As with NS/NS, the coalescence 
waves should carry equation of state information and will come out in the kHz 
band, where their detection will require advanced, specialty detectors. 



Christodoulou Memory 

As the coalescence waves depart from their source, their energy creates (via 
the nonlinearity of Einstein's field equations) a secondary wave called the "Chris- 



todoulou memory" |5^, |5^, Whereas the primary waves may have frequencies 
in the kHz band, the memory builds up on the timescale of the primary energy 
emission profile, which is likely to be of order 0.01 sec, corresponding to a memory 
frequency in the optimal band for the LIGO/VIRGO workhorse interferometers, 
~ lOOHz. Unfortunately, the memory is so weak that only very advanced ground- 
based interferometers have much chance of detecting and studying it — and then, 
perhaps only for BH/BH coalescences and not for NS/NS or NS/BH [56|. LISA, by 



contrast, should easily be able to measure the memory from supermassive BH/BH 
coalescences. 



5 Other High- Frequency Sources 
5.1 Stellar Core Collapse and Supernovae 

When the core of a massive star has exhausted its supply of nuclear fuel, it col- 
lapses to form a neutron star or black hole. In some cases, the collapse triggers 
and powers a subsequent explosion of the star's mantle — a supernova explosion. 
Despite extensive theoretical efforts for more than 30 years, and despite wonderful 
observational data from Supernova 1987A, theorists are still far from a definitive 
understanding of the details of the collapse and explosion. The details are highly 
complex and may differ greatly from one core collapse to another |57| |. 

Several features of the collapse and the core's subsequent evolution can produce 
significant gravitational radiation in the high-frequency band. We shall consider 
these features in turn, the most weakly radiating first. 

Boiling of the Newborn Neutron Star 

Even if the collapse is spherical, so it cannot radiate any gravitational waves at 
all, it should produce a convectively unstable neutron star that "boils" vigorously 
(and nonspherically) for the first ~ 1 second of its life The boiling dredges 

up high-temperature nuclear matter (T ~ lO^^K) from the neutron star's central 
regions, bringing it to the surface (to the "neutrino-sphere"), where it cools by 
neutrino emission before being swept back downward and reheated. Burrows |l59|| 
has pointed out that the boiling should generate n ~ 100 cycles of gravitational 
waves with frequency / ~ lOOHz and amplitude large enough to be detectable 
by LIGO/VIRGO throughout our galaxy and its satellites. Neutrino detectors 
have a similar range, and there could be a high scientific payoff from correlated 
observations of the gravitational waves emitted by the boiling's mass motions and 
neutrinos emitted from the boiling neutrino-sphere. With neutrinos to trigger on, 



the sensitivities of LIGO detectors should be about twice as good as shown in Fig. 
I 

Recent 3+1 dimensional simulations by Miiller and Janka |6C] suggest an rms 



amplitude /i ~ 2 x 10 ^^(20kpc/r) (where r is the distance to the source), corre- 
sponding to a characteristic amplitude he — hy/n ~ 2 X 10-22 (20kpc/r); cf. Fig. 

(The older 2+1 dimensional simulations give he about 6 times larger than this 
[]60| , but presumably are less reliable.) LIGO should be able to detect such waves 
throughout our galaxy with an amplitude signal to noise ratio of about S/N = 2.5 
in each of its two enhanced 4km interferometers, and its advanced interferometers 
should do the same out to 80Mpc distance. (Recall that the /isB curves in Fig. § are 
drawn at a signal to noise ratio of about 5). Although the event rate is only about 
one every 40 years in our galaxy and one every 20 years out to 80Mpc, the cor- 
related neutrino and gravitational wave observations could bring very interesting 
insights into the boiling of a newborn neutron star. 

Axisymmetric Collapse, Bounce, and Oscillations 

Rotation will centrifugally flatten the collapsing core, enabling it to radiate 
as it implodes. If the core's angular momentum is small enough that centrifugal 
forces do not halt or strongly slow the collapse before it reaches nuclear densities, 
then the core's collapse, bounce, and subsequent oscillations are likely to be axially 
symmetric. Numerical simulations |61, |63] show that in this case the waves from 



collapse, bounce, and oscillation will be quite weak: the total energy radiated as 
gravitational waves is not likely to exceed ~ 10"^ solar masses (about 1 part in a 
million of the collapse energy) and might often be much less than this; and corre- 
spondingly, the waves' characteristic amplitude will he he ^ 3 x 10~2i(30kpc/r). 
These collapse-and-bounce waves will come off at frequencies ~ 200 Hz to ~ 1000 
Hz, and will precede the boiling waves by a fraction of a second. Though a little 
stronger than the boiling waves, they probably cannot be seen by LIGO/VIRGO 
beyond the local group of galaxies and thus will be a very rare occurrence. 



Rotation- Induced Bars and Break- Up 

If the core's rotation is large enough to strongly flatten the core before or as 
it reaches nuclear density, then a dynamical and/or secular instability is likely to 
break the core's axisymmetry. The core will be transformed into a bar-like con- 
figuration that spins end-over-end like an American football, and that might even 
break up into two or more massive pieces. In this case, the radiation from the 
spinning bar or orbiting pieces could be almost as strong as that from a coalescing 
neutron-star binary, and thus could be seen by the LIGO/VIRGO first interferom- 
eters out to the distance of the Virgo cluster (where the supernova rate is several 
per year), by enhanced interferometers out to ~ lOOMpc (supernova rate several 
thousand per year), and by advanced interferometers out to several hundred Mpc 
(supernova rate ~ (a few)x ~ 10^ per year); cf. Fig. It is far from clear what 
fraction of collapsing cores will have enough angular momentum to break their ax- 
isymmetry, and what fraction of those will actually radiate at this high rate; but 
even if only ~ 1/1000 or 1/10^ do so, this could ultimately be a very interesting 
source for LIGO/VIRGO. 



Several specific scenarios for such non-axisymmetry have been identified: 
Centrifugal hangup at ~ 100km radius: If the pre-collapse core is rapidly 
spinning (e.g., if it is a white dwarf that has been spun up by accretion from a com- 
panion), then the collapse may produce a highly flattened, centrifugally supported 
disk with most of its mass at radii R ~ 100km, which then (via instability) may 
transform itself into a bar or may bifurcate. The bar or bifurcated lumps will radi- 
ate gravitational waves at twice their rotation frequency, / ~ lOOHz — the optimal 
frequency for LIGO/VIRGO interferometers. To shrink on down to ~ 10km size, 
this configuration must shed most of its angular momentum. If a substantial frac- 
tion of the angular momentum goes into gravitational waves, then independently 
of the strength of the bar, the waves will be nearly as strong as those from a coa- 
lescing binary. The reason is this: The waves' amplitude h is proportional to the 
bar's ellipticity e, the number of cycles n of wave emission is proportional to 1/e^, 
and the characteristic amplitude he = h^/n is thus independent of the ellipticity 
and is about the same whether the configuration is a bar or is two lumps |31]. The 
resulting waves will thus have he roughly half as large, at / ~ lOOHz, as the he 
from a NS/NS binary (half as large because each lump might be half as massive as 
a NS), and the waves will chirp upward in frequency in a manner similar to those 
from a binary. 

It is rather likely, however, that most of excess angular momentum does not go 
into gravitational waves, but instead goes largely into hydrodynamic waves as the 
bar or lumps, acting like a propeller, stir up the surrounding stellar mantle. In this 
case, the radiation will be correspondingly weaker. 

Centrifugal hangup at ~ 20km radius: Lai and Shapiro [^3| have explored 
the case of centrifugal hangup at radii not much larger than the final neutron star, 
say R ~ 20km. Using compressible ellipsoidal models, they have deduced that, after 
a brief period of dynamical bar-mode instability with wave emission at / ~ lOOOHz 
(explored by Houser, Centrella, and Smith |64]), the star switches to a secular 
instability in which the bar's angular velocity gradually slows while the material of 
which it is made retains its high rotation speed and circulates through the slowing 
bar. The slowing bar emits waves that sweep downward in frequency through 
the LIGO/VIRGO optimal band / ^ lOOHz, toward ~ lOHz. The characteristic 
amplitude (Fig. ^ is only modestly smaller than for the upward-sweeping waves 
from hangup at i? ~ lOOkm, and thus such waves should be detectable near the 
Virgo Cluster by the first LIGO/VIRGO interferometers, near lOOMpc by enhanced 
interferometers, and and at distances of a few lOOMpc by advanced interferometers. 

Successive fragmentations of an accreting, newborn neutron star: 
Bonnell and Pringle ||65[| have focused on the evolution of the rapidly spinning, 
newborn neutron star as it quickly accretes more and more mass from the pre- 
supernova star's inner mantle. If the accreting material carries high angular mo- 
mentum, it may trigger a renewed bar formation, lump formation, wave emission, 
and coalescence, followed by more accretion, bar and lump formation, wave emis- 
sion, and coalescence. Bonnell and Pringle speculate that hydrodynamics, not wave 
emission, will drive this evolution, but that the total energy going into gravitational 
waves might be as large as ~ 1O~^M0. This corresponds to he ~ 10~^^(10Mpc/r). 



5.2 Spinning Neutron Stars; Pulsars 



As the neutron star settles down into its final state, its crust begins to solidify 
(crystalize). The solid crust will assume nearly the oblate axisymmetric shape 
that centrifugal forces are trying to maintain, with poloidal ellipticity ep oc (angular 
velocity of rotation)^. However, the principal axis of the star's moment of inertia 
tensor may deviate from its spin axis by some small "wobble angle" Oyj, and the 
star may deviate slightly from axisymmetry about its principal axis; i.e., it may 
have a slight ellipticity ee <C in its equatorial plane. 



As this slightly imperfect crust spins, it will radiate gravitational waves |66|: 
ee radiates at twice the rotation frequency, / = 2/i.ot with h oc €e, and the wobble 
angle couples to ep to produce waves at / = /rot+/prec (the precessional sideband of 
the rotation frequency) with amplitude h oc Oyjdp. For typical neutron-star masses 
and moments of inertia, the wave amplitudes are 



[^) i'^) . (3) 



The neutron star gradually spins down, due in part to gravitational-wave emis- 
sion but perhaps more strongly due to electromagnetic torques associated with its 
spinning magnetic field and pulsar emission. This spin-down reduces the strength 
of centrifugal forces, and thereby causes the star's poloidal ellipticity ep to decrease, 
with an accompanying breakage and resolidification of its crust's crystal structure 
(a "starquake") [^]. In each starquake, 6^, ee, and ep will all change suddenly, 
thereby changing the amplitudes and frequencies of the star's two gravitational 
"spectral lines" / = 2/rot and / = /rot + /prec- After each quake, there should 
be a healing period in which the star's fluid core and solid crust, now rotating 
at different speeds, gradually regain synchronism. By monitoring the amplitudes, 
frequencies, and phases of the two gravitational-wave spectral lines, and by com- 
paring with timing of the electromagnetic pulsar emission, one might learn much 
about the physics of the neutron-star interior. 

How large will the quantities ee and Ow^p be? Rough estimates of the crustal 
shear moduli and breaking strengths suggest an upper limit in the range e^ax ^ 
10~^ to 10~^, and it might be that typical values are far below this. We are 
extremely ignorant, and correspondingly there is much to be learned from searches 
for gravitational waves from spinning neutron stars. 

One can estimate the sensitivity of LIGO/VIRGO (or any other broad-band 
detector) to the periodic waves from such a source by multiplying the waves' am- 
plitude h by the square root of the number of cycles over which one might integrate 
to find the signal, n = ff where f is the integration time. The resulting effective 
signal strength, h^/n, is larger than h by 



lOOOHz/ V4monthsy 

Four months of integration is not unreasonable in targeted searches; but for an 
all-sky, all-frequency search, a coherent integration might not last longer than a 
few days because of computational limitations associated with having to apply 



huge numbers of trial neutron-star spindown corrections and earth-motion doppler 
corrections |]68[| . 

Equation (Q) for h^/n should be compared (i) to the detector's rms broad-band 
noise level for sources in a random direction, \/5^rms! to deduce a signal-to-noise 
ratio, or (ii) to hsB to deduce a sensitivity for high-confidence detection when one 
does not know the waves' frequency in advance Q. Such a comparison suggests that 
the first interferometers in LIGO/VIRGO might possibly see waves from nearby 
spinning neutron stars, but the odds of success are very unclear. 

The deepest searches for these nearly periodic waves will be performed by 
narrow-band detectors, whose sensitivities are enhanced near some chosen fre- 
quency at the price of sensitivity loss elsewhere — e.g., dual-recycled interferometers 



[15 1 or resonant-mass antennas (Section 2A). With "advanced-detector technology" 



and targeted searches, dual-recycled interferometers might be able to detect with 
confidence spinning neutron stars that have ||^] 

(eeor0^ep)^3xio-i° — tt^ttt^ ■ (5) 



V /rot J VlOOOpc 

There may well be a large number of such neutron stars in our galaxy; but it is 
also conceivable that there are none. We are extremely ignorant. 

Some cause for optimism arises from several physical mechanisms that might 
generate radiating ellipticities large compared to 3 x 10"^'': 

• It may be that, inside the superconducting cores of many neutron stars, 
there are trapped magnetic fields with mean strength Score ~ 10^^ G or even 
lO^^G. Because such a field is actually concentrated in flux tubes with B = 
-Bcrit ~ 6 X 10^^ G surrounded by fleld-free superconductor, its mean pressure 
is pb = -Bcorc-Bcrit/Svr. This pressure could produce a radiating ellipticity eg ~ 
^toCp ^ Pb/p ^ lO^^-Bcorc/lO^'^^G (where p is the core's material pressure). 

• Accretion onto a spinning neutron star can drive precession (keeping 9^] sub- 



stantially nonzero), and thereby might produce measurably strong waves [69|. 



If a neutron star is born rotating very rapidly, then it may experience a 
gravitational-radiation-reaction-driven instability. In this "CFS" (Chandra- 
sekhar, |7C] Friedman, Schutz [^]) instability, density waves travel around 



the star in the opposite direction to its rotation, but are dragged forward 
by the rotation. These density waves produce gravitational waves that carry 
positive energy as seen by observers far from the star, but negative energy 
from the star's viewpoint; and because the star thinks it is losing negative 
energy, its density waves get amplifled. This intriguing mechanism is similar 
to that by which spiral density waves are produced in galaxies. Although the 
CFS instability was once thought ubiquitous for spinning stars f% we 



now know that neutron-star viscosity will kill it, stabilizing the star and turn- 
ing off the waves, when the star's temperature is above some limit ~ lO^'^K 
|73| and below some limit ~ lO^K |74|; and correspondingly, the instability 
should operate only during the first few years of a neutron star's life, when 
lO^K <,T<, IQiOR. 



6 Low-Frequency Gravitational- Wave Sources 



6.1 Waves from the Coalescence of Massive Black Holes in Distant 
Galaxies 

LISA would be a powerful instrument for studying massive black holes in distant 
galaxies. Figure ^ shows, as examples, the waves from several massive black hole 
binaries at 3Gpc distance from Earth (a cosmological redshift of unity). The waves 
sweep upward in frequency (rightward in the diagram) as the holes spiral together. 
The black dots show the waves' frequency one year before the holes' final collision 
and coalescence, and the arrowed lines show the sweep of frequency and char- 
acteristic amplitude he = h^/n during that last year. For simplicity, the figure is 
restricted to binaries with equal-mass black holes: W'^Mq/W'^Mq, IO^Mq/IO^Mq, 
and W^Mq/W^Mq. 

By extrapolation from these three examples, we see that LISA can study much 
of the last year of inspiral, and the waves from the final collision and coalescence, 
whenever the holes' masses are in the range 3 x W^Mq ^ M ^ 3 x W^Mq p6| . 
Moreover, LISA can study the final coalescences with remarkable signal to noise 
ratios: S/N ^ 1000. Since these are much larger 5/iV's than LIGO/VIRGO is 
likely to achieve, we can expect LISA to refine the experimental understanding of 
black-hole physics, and of highly nonlinear vibrations of warped spacetime, which 
LIGO /VIRGO initiates — provided the rate of massive black- hole coalescences is of 
order one per year in the Universe or higher. The rate might well be that high, 
but it also might be much lower. 

By extrapolating Fig. |6| to lower BH/BH masses, we see that LISA can observe 
the last few years of inspiral, but not the final collisions, of binary black holes in 



the range lOOM© <,M ^ 10^ Mq, out to cosmological distances 

Extrapolating the BH/BH curves to lower frequencies using the formula (time to 
final coalescence) oc f^^^^, we see that equal-mass BH/BH binaries enter LISA's 
frequency band roughly 1000 years before their final coalescences, more or less 
independently of their masses, for the range IOOMq <^ M <, 10^ Mq. Thus, if 
the coalescence rate were to turn out to be one per year, LISA would see roughly 
1000 additional massive binaries that are slowly spiraling inward, with inspiral 
rates df /dt readily measurable. From the inspiral rates, the amplitudes of the two 
polarizations, and the waves' harmonic content, LISA can determine each such 
binary's luminosity distance, redshifted chirp mass (1 + 2;)Mc, orbital inclination, 
and eccentricity; and from the waves' modulation by LISA's orbital motion, LISA 
can learn the direction to the binary with an accuracy of order one degree. 



6.2 Waves from Compact Bodies Spiraling into Massive Black 
Holes or Exotic Objects in Distant Galaxies 

When a compact body with mass ^ spirals into a much more massive black hole with 
mass M, the body's orbital energy E at fixed frequency / (and correspondingly at 
fixed orbital radius a) scales as E oc fi, the gravitational- wave luminosity E scales 
as E on y?, and the time to final coalescence thus scales as t ~ E/E oc l//i. This 
means that the smaller is /i/M, the more orbits are spent in the hole's strong- 
gravity region, a <, 10GM/c^ and thus the more detailed and accurate will be the 



map of the hole's spacetime geometry, which is encoded in the emitted waves. 

For holes observed by LIGO /VIRGO, the most extreme mass ratio that we 
can hope for is ^i/M ~ IMq/SOOM©, since for M > SOOMq the inspiral waves 
are pushed to frequencies below the LIGO/VIRGO band. This limit on n/M 
seriously constrains the accuracy with which LIGO/VIRGO can hope to map out 
the spacetime geometries of black holes and test the black-hole no-hair theorem 
[48, 49 1 (end of Section ^^). By contrast, LISA can observe the final inspiral waves 
from objects of any mass M ^ 0.5Mq spiraling into holes of mass 3 x 10^ Mq ^ 
M ^ 3 X IO^Mq. 

Figure |6| shows the example of a IOMq black hole spiraling into a 10^ Mq hole at 
3Gpc distance. The inspiral orbit and waves are strongly influenced by the hole's 
spin. Two cases are shown [75|: an inspiraling circular orbit around a non-spinning 
hole, and a prograde, circular, equatorial orbit around a maximally spinning hole. 
In each case the dot at the upper left end of the arrowed curve is the frequency and 
characteristic amplitude one year before the final coalescence. In the nonspinning 
case, the small hole spends its last year spiraling inward from r ~ 7AGM/c^ (3.7 
Schwarzschild radii) to its last stable circular orbit at r = GGM/c^ (3 Schwarzschild 
radii). In the maximal spin case, the last year is spent traveling from r = QGM/(? 
(3 Schwarzschild radii) to the last stable orbit at r = GM/c? (half a Schwarzschild 
radius). The ~ 10^ cycles of waves during this last year should carry, encoded in 
themselves, rather accurate values for the massive hole's lowest few multipole mo- 
ments ||48[| (or, equivalently, a fairly accurate map of the hole's spacetime geometry 
(or, hole's). If the measured moments satisfy the "no-hair" theorem (i.e., if they are 
all determined uniquely by the measured mass and spin in the manner of the Kerr 
metric), then we can be sure the central body is a black hole. If they violate the 
no-hair theorem, then (assuming general relativity is correct), either the central 
body was an exotic object (e.g. soliton star or naked singularity) rather than a 
black hole, or else an accretion disk or other material was perturbing its orbit |76|. 
From the evolution of the waves one can hope to determine which is the case, and 
to explore the properties of the central body and its environment |p9t| . 

Models of galactic nuclei, where massive holes reside, suggest that inspiraling 
stars and small holes typically will be in rather eccentric orbits [^]. This is because 
they get injected into such orbits via gravitational deflections off other stars, and 
by the time gravitational radiation reaction becomes the dominant orbital driving 
force, there is not enough inspiral left to fully circularize their orbits. Such or- 
bital eccentricity will complicate the waveforms and complicate the extraction of 
information from them. Efforts to understand the emitted waveforms are just now 
getting underway. 

The event rates for inspiral into massive black holes are not at all well under- 
stood. However, since a significant fraction of all galactic nuclei are thought to 
contain massive holes, and since white dwarfs and neutron stars, as well as small 
black holes, can withstand tidal disruption as they plunge toward the massive hole's 
horizon, and since LISA can see inspiraling bodies as small as ~ O.SM© out to 3Gpc 
distance, the event rate is likely to be interestingly large. 



7 Conclusion 



It is now 36 years since Joseph Weber initiated his pioneering development of 



gravitational- wave detectors y] and 25 years since Forward and Weiss |79| ini- 
tiated work on interferometric detectors. Since then, hundreds of talented experi- 
mental physicists have struggled to improve the sensitivities of these instruments. 
At last, success is in sight. If the source estimates described in this lecture are ap- 
proximately correct, then the planned interferometers should detect the first waves 
in 2001 or several years thereafter, thereby opening up this rich new window onto 
the Universe. 
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